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' Abstract 

QQ ' The Global Positioning System (GPS) provides an excellent educa- 

tional example as to how the theory of general relativity is put into prac- 
■ tice and becomes part of our everyday life. This paper gives a short and 

^ ' instructive derivation of an important formula used in the GPS, and is 

aimed at graduate students and general physicists. 
QP ■ The theoretical background of the GPS (see [1]) uses the Schwarzschild 

spacetime to deduce the approximate formula , ds/dt ~ 1 + V — , for the 
relation between the proper time rate s of a satellite clock and the coordi- 
' nate time rate t. Here V is the gravitational potential at the position of the 

, satellite and v is its velocity (with light-speed being normalized as c = 1). 

In this note we give a different derivation of this formula, without us- 



^ ^, ing approximations, to arrive at ds/dt = y^l -f 2V — |vp — ^^^^ (n ■ v)^, 

where n is the normal vector pointing outward from the center of Earth 
to the satellite. In particular, if the satellite moves along a circular orbit 
' then the formula simplifies to ds/dt = \/l + 2V — |vp. 

^ I We emphasize that this derivation is useful mainly for educational 

• • . purposes, as the approximation above is already satisfactory in practice. 

> 



1 Introduction 

The most significant application of the theory of General Relativity in everyday 
life, arguably, is the Global Positioning System. The GPS uses accurate, stable 
atomic clocks in satellites and on the ground to provide world-wide position 
and time determination. These clocks have relativistic frequency shifts which 
need to be carefully accounted for, in order to achieve synchronization in an 
underlying Earth-centered inertial frame, upon which the whole system is based. 
For educational purposes it is very instructive to see how a highly abstract theory 
such as the general theory of relativity becomes a part of our everyday life. 
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Throughout this note we will normalize universal constants for simplicity, so 
that light speed and the gravitational constant are the unity, c = 1, 7 = 1. 

Let us briefly sketch here how the GPS works, including the major ingredi- 
ents for the reader's convenience (a full and detailed description is available in 
[T]). In order to determine your position at the surface of Earth you must be 
able to see some (at least four) GPS-satellites simultaneously, and know their 
position and your distance to them (with this much information at hand, it is 
then elementary geometry to determine your position). To make this possible, 
every GPS-satellite emits signals continuously, describing its position and its 
local time tsat- If your GPS-device measures its own local time, tdev, then you 
can use the constancy of the speed of light c to calculate your distance d to the 
satellite as 

^ {tdev isat)c — ^dev tsat (1) 

keeping in mind the normalization c = 1. (We remark here that in practice the 
GPS-device is unable to measure time with sufficient accuracy, therefore tdev 
also has to be deduced from the data sent by the satellites. However, this detail 
is not important for the purposes of this paper.) 

This looks simple enough, but the problem is that formula ([T]) is only valid if 
all clocks are synchronized in some special relativistic inertial frame. Therefore, 
for the purposes of the GPS one imagines that an inertial frame is attached to 
the center of Earth, and we try to synchronize all clocks such that they measure 
the time t of this ideal inertial frame. This means that one should imagine 
'ideal' clocks placed everywhere in the vicinity of Earth, measuring the time t, 
and thus tsat and tdev should be the read-outs of these ideal clocks at the place 
of the satellite when the signal originates and at the place of your device when 
the signal is received. However, what we actually can measure is the proper 
time tE oi stationary clocks on the surface of Earth (which is measured in time- 
keeping centers throughout the world), and the proper time s of satellite clocks. 
Therefore, we need to establish a relation between the time rates tE, s and t. 

Fortunately, the time rate Ie measured on Earth is independent of where you 
are (i.e. the same manufactured clocks beat the same rate in London or Tokyo 
or New York) and differs from the 'ideal' time rate t only by a multiplicative 
constant, 

;^ = l-<i>o, (2) 
dtE 

where $0 is a constant corresponding to the Earth's geoid (see [1]). This re- 
lation is very convenient because the ideal time rate t can be replaced by t^, 
something we can actually measure, and an equation of the form ([Ij still re- 
mains valid after rescaling by the factor 1 — $o- This leaves us with the task of 
establishing a relation between s and t (or s and tE, whichever turns out to be 
more convenient). 

To determine the relation of the proper time rate s of the satellite clocks and 
the time rate t of ideal clocks measuring the coordinate-time of the underlying 
Earth-centered inertial frame, the customary theoretical framework [IJ is to use 
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Schwarzschild spacetime, and arrive at the formula 

IvP 

ds/dt Kl + V - (3) 

after several first-order approximations in the calculations. Here V denotes the 
gravitational potential at the position of the satellite, and v is the velocity of 
the satellite measured in the underlying non-rotating Earth-centered inertial 
frame. Formula ^ is the internationally accepted standard relating the clock 
frequencies, as described in [H [3] and references therein. 

We remark here that the derivation of formula ^ is somewhat more involved 
than that of formula ([3]). In fact, for the purposes of deriving ^ one needs to 
modify the Schwarzschild metric by a small term, taking into account multipole 
contributions corresponding to the Earth's geoid. However, for formula ([3]) the 
modifying term is disregarded and the standard Schwarzschild metric is used, 
the argument being that the modifying term becomes negligible far enough from 
the Earth surface, where the satellites orbit (cf. [1]). In this paper we do not 
include the derivation of formula ([2]) and the modifying term which is used there. 

In general, in physics it is justified to use approximate formulae for two 
different reasons. One is that in some cases the derivation of an exact result is 
analytically not possible. The other is that in certain cases an exact analytic 
derivation, even if it exists, would lead to involved and lengthy calculations thus 
concealing the important and possibly very simple aspects at the heart of the 
issue. In calculations involving the theory of relativity (in particular, general 
relativity) there is a tendency to turn to approximations automatically, due to 
the involved nature of the theory. However, in some rare cases a modified point 
of view and an adequate choice of coordinate system can lead to exact results. 

In this note, we adhere to the standard theoretical framework used in [TJ[3], 
but we point out that an exact formula can be derived in a very simple manner 
for the clock frequency rate ^ in question. Instead of using the customary 
isotropic coordinates we will treat Schwarzschild spacetime with Schwarzschild 
coordinates (an entirely coordinate-free derivation of the same formulae can be 
found in [5^ but it is rather cumbersome). This treatment of Schwarzschild 
spacetime is motivated by a similar account of special relativity in ^ and that 
of general relativity in 7 . The coordinate-free point of view often has the ad- 
vantage of conceptual clarity and, in this particular case, brevity of calculations. 

Our results here are mostly of educational and theoretical interest as the 
existing formula ([3]) provides good approximation to the desired precision in 
the GPS (see [3]). 

2 Schwarzschild 's spacetime 

Schwarzschild's spacetime describes the gravitational field of a pointlike inertial 
mass TO. It is a well-known model of general relativity, but we include its short 
description here for convenience. 
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Let us introduce some notation. Let E he a. three-dimensional Euclidean 
vector space, the inner product of x, y S E being denoted by x-y. For x € E 
we put 

n(x) := ^ (4) 

for the outward normal vector at x. 

Consider R x as a spacetime manifold with its usual special relativistic 
metric, i.e. the Lorentz form given by 

where / is the identity matrix. For a world point (t, x) E R x E one should 
think of t as the synchronization time corresponding to the center of Earth, and 
X as the space vector pointing to the world point from the center of Earth. The 
Lorentz form (O means that the Lorentz length-square of a four-vector (s,q) 
is given by the usual formula |(s,q)||Q,. = — + |qp. This special relativistic 
spacetime will be called the Earth Centered Reference Frame (ECRF). This 
is an 'ideal' special relativistic frame, and the task in the GPS is to achieve 
synchronization of clocks in this underlying inertial frame. 

Introduce the potential V{x) — on i?, and restrict your consideration 
to world points (t, x) E R x E for which 1 -I- 2V{x) > (i.e. for world points 
outside the Schwarzschild radius). Then it is easy to see (cf. [SIEIIH]) that for 
such world points, the standard form of the Schwarzschild metric in Rx _B (i.e. a 
smooth collection of Lorentz forms g{t, x) depending on the world points) takes 
the form: 



,-(l + 2T/(x)) , 

^^''^^^=1 /-^n(x)«n(x).J 



In particular, when we measure the length-square of any space vector q at 
the space point x in Schwarzschild's metric we obtain 



NL.-|qp-^f^(n(x).q))l (7) 

Similarly, when we measure the Schwarzschild length-square of a four-vector 
(e.g. a four-velocity) (s,q) at the point x we get 

l(s,q)lL. = -(i + 2nx))).' + |q|L.. (8) 



3 A satellite in Earth's gravitational field 

In this section we derive the formula relating the proper time s of satellite 
vehicle-clocks to the ideal time t of the underlying inertial frame. 

A material point in spacetime is described by a world line function 

R ^ M X s {t{s),x{s)) (9) 



4 



where s is the proper time of the material point. The four-velocity, (i(s), i{s)) , 
of the material point always satisfies 

\{i{s),x{s))\l,^ = -l. (10) 

Therefore, using equation ([5]) we obtain 

- (1 + 2V{x{s)))i{sf + |i(s)||,;, = -1 (11) 

The time function t{s) can be inverted, with s{t) denoting the proper time 
instant s corresponding to the ECRF-time t. Let v{t) := ^^^^ denote the 
relative velocity of the material point with respect to the center of the Earth. 
Note that by the chain rule we have v(t) = x{s{t))j^. Now, rearranging (jlip 
we obtain 

(f)%ii(.(t))iL,(i)!^i (12) 



l + 2Vix{s{t))) 
which yields 



^ = ^l + 2V{x{sm-\v{t)\l^, (13) 
Finally, applying ^ we obtain the desired relation 



— = \1 + 2V-M^ ?^(n.v)2 (14) 

Using a series expansion, assuming that all the terms on the right-hand side 
are much less than 1, we get back the approximate formula ^ 

Of course, there is still some work to be done before one can apply equation 
^ in the GPS. Instead of ^ what we really need is the function t{s), because 
the signal emitted by the satellite contains the proper time instant s measured 
by clocks on board, and we would like to replace it by the coordinate-time 
instant t which it corresponds to. However, having deduced formula (|14[) (or 
its approximation Q) it is possible to obtain the function i(s), or at least a 
good enough approximation of it. This is described very well in full detail 
in [1], and we feel it inappropriate to repeat those calculations word- by- word 
here. Nevertheless, we warmly recommend that the reader turn to ,lj for the 
interesting details. 
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